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ABSTRACT 

The unconditional mass function (UMF) of dark matter haloes has been determined 
accurately in the literature, showing excellent agreement with high resolution numeri- 
cal simulations. However, this is not the case for the conditional mass function (CMF). 
Here, we propose a simple analytical procedure to derive the CMF by rescaling the 
UMF to the constrained environment using the appropriate mean and variance of the 
density field at the constrained point. This method introduces two major modifications 
with respect to the standard re-scaling procedure. First of all, rather than using in the 
scaling procedure the properties of the environment averaged over all the conditioning 
region, we implement the re-scaling locally. We show that for high masses this mod- 
ification may lead to substantially different results. Secondly, we modify the (local) 
standard re-scaling procedure in such a manner as to force normalisation, in the sense 
that when one integrates the CMF over all possible values of the constraint multiplied 
by their corresponding probability distribution, the UMF is recovered. In practise, we 
do this by replacing in the standard procedure the value <5 C (the linear density contrast 
for collapse) by certain adjustable effective parameter S e s. In order to test the method, 
we co mpare our prescr i ption with the results obtained from numerical simulations in 
voids (|Gottl5ber et al.ll2003D . We find that when our modified re-scaling is applied 
locally to any existing numerical fit of the UMF, and the appropriate value for (5 e ff is 
chosen, the resulting CMF is, in all cases, in very good agreement with the numerical 
results. Based on these results, we finally present a very accurate analytical fit to the 
(accumulated) conditional mass function obtained with our procedure, as a function 
of the parameters that describe the conditioning region (size and mean linear density 
contrast), the redshift and the relevant cosmological parameters (as and T). This ana- 
lytical fit may be useful for any theoretical treatment of the large scale structure, and 
has been already used successfully in regard with the statistic of voids. 



Key words: methods: analytical - methods: statistical 
matter - large-scale structure of the Universe 



cosmology: theory - dark 



1 INTRODUCTION 

In the last few years, there has been considerable effort 
in obtaining accurate theoretical predictions for the mass 
function of collapsed dark matter haloes. By far, the most 
widely used predi ction for the unco n dition al mass func- 
tion (UMF) is the IPress fc Schechterl l|l974 hereafter PS) 
formalism. Extensions of this prescription provide a way 
to compute not only good approximation s to the UMF 
jPress fc Schechterl fl97l IBond et al l Il99ll ). but also the 
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merging history (IBond et al. 199ll: Lacev fc Cold 199&t) and 
the s patial clustering l|Mo fc White! 1 19961 : ISheth fc Lemsonl 
1999) of dark matter haloes. 

The two basic assumptions of the PS approach are the 
physics of the spherical collapse, and the fact that the ini- 
tial fluctua ti ons w ere drawn from a gaussian distribution. 
I Bond et al.l l|l99ll ) showed how to combine these two hy- 
pothesis in order to obtain the UMF of dark matter haloes 
from the barrier crossing statistics of many independent, 
uncorrelated random walks (the so-called "excursion set for- 
malism"). The barrier shape is given by the fact that, in the 
spherical collapse model, a certain region collapses at time 
t (or redshift z) if the initial overdensity within it exceeds 
a critical value (5 SC ) which is independent of mass. How- 
ever, numerical simulations show that the PS UMF as de- 
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rived from this formalism, while qualitatively correct, over- 
estimates the abundance of "typical" haloes and underesti- 
mates that of massive ones when compared with t he results 
of N-body simulations (e.g. lEfstathiou et al.|[l988h . 

In order to understand if this discrepancy could be due 
to the assumption of spherical collapse made within the 
PS formalism, several works in the last few years extended 
the excursion set formali sm by incorporating a treatment 
of th e ellipsoidal collapse (jBond fc Mverdll996l : ISheth et al] 
l200ll : ISheth fc Tormenl I2OO2I ). In practise, this is done by 
using a "moving barrier", i.e. a barrier whose height de- 
pends on mass. This study provides an analytic expression 
which describes to a good approximation the distribution 
of first-crossings of this "moving" barrier, once two free pa- 
rameters are fitted by c omparison with N-body simulations 
|Sheth fc Tormenl Il999l . ST). The resulting mass function 
significantly improves the results of the PS formalism, al- 
though a small disc repancy still remains at high masses 
l|jenkins et al-l feoOl), The authors of this last paper also 
propose an analytic fitting formula which better reproduces 
the numerical results, but which can not be extrapolated in 
mass beyond the range of the fit. 

Recently, IWarren et all l|2006l . WA) use a large sample 
of simulations to provide a extrem ely good fitt i ng fo rmula 
for the UMF at redshift zero, while iReed et ail (120071 ) pro- 
vide a good fit to the UMF at redshifts 1 0-30. O n the other 
hand. iBetancort-Riio fc Montero-Dorta ( 2006a. b) proposed 
an analytical procedure for deriving the UMF at any red- 
shift, which also includes the physics of the ellipsoidal col- 
lapse encoded in the ST mass function. In this procedure, 
the "all-mass-at-center" problem is treated in an appropri- 
ate manner. As a consequence of this, the high mass be- 
haviour and the redshift dependence of the UMF comes out 
right without introducing any fudge parameter to be fit- 
ted with simulation s. Their expression reproduces well the 
IWarren et ail (2006) fit to the simulations. 

Summarising, there is a variety of analytical expressions 
which accurately reproduce the UMF. However, this is not 
the case for the conditional mass function (CMF). There are 
different approaches in the literature to obtain an analytic 
expression for the CMF of dark matter haloes, n c (m). The 
most widely used framework is again the excur sion set for- 
malism, often called exte nded Press-Schechter (|Bond et al.l 
Il99ll ; iLacey fc Cole Il993l ). and denoted EPS (for a recent 
review, see Izentnedl2007f ). However, this appro ach does no t 
compares well to N-body simulations (see e.g. iTormenll 199SI ) . 
being one of the reasons t hat it does not include th e physics 
of the ellipsoidal collapse. ISheth fc Tormenl (|2002T l propose 
an alternative expression of the CMF to solve this issue 
within the context of the ellipsoidal collapse moving barrier, 
which reproduces much better the results from simulations 
of the hierarchical assembly of dark matter haloes. 

In this paper, we present our method to build the CMF 
of dark matter haloes. Our motivation is to provide an an- 
alytical prescription to build the CMF as an extension of 
the UMF, taking advantage of the high degree of accuracy 
at which we know these functions. Our proposal consists in 
a prescription to rescale the UMF, which can be applied 
to any of the existing fitting formula, and which takes into 
account a certain normalisation condition that we will de- 
scribe below. Our procedure differs from the standard one in 
two major respects. First, the scaling of the UMF is imple- 



mented locally, using at each point within the conditioning 
region the local environment, rather than using the average 
(in practise, this means that our CMF will depend on an 
additional variable, q, the distance from the center of the 
condition to the point at which we evaluate the CMF, and 
over which we have to integrate in order to obtain the av- 
erage CMF). Secondly, the scaling procedure is somewhat 
different from the standard one, and instead of the usual 
value of 8 C , we will use a modified one. 

We assume throughout the paper a flat ACDM model 
with the following cosmological parameters: f2 m = 0.3, 
h = 0.7 (i.e. F = n m h = 0.21), cr 8 = 0.9 and no tilt of 
the primordial spectrum (i.e. a scalar spectral index ns of 
1). The CDM power spectrum (P(k) oc k ns T(k) 2 ) adopted 
for our compu tations is obtained f r om th e transfer function 
T(k) given bv lBond fc Efstathioul i| 19841 ). 



T(k) = ( 1 + [ak + {bkf /2 + (ck) 2 }' 



-l/u 



(1) 



where a = 6.4/F h" 1 Mpc, b 
1.7/r h" 1 Mpc, and v = 1.13. 



3.0/F h _1 Mpc, c 



2 THE UNCONDITIONAL MASS FUNCTION 

The (unconditional) mass function (UMF) of dark matter 
haloes at a given redshift z, n(m,z), is defined such that 
n(m, z)dm is the comoving number density of haloes in the 
mass range (m, m + dm) at that redshift. The bound ob- 
jects (haloes) are usually defined in numerical simulations 
using two di f ferent algorithms: the friends-of-friends finder 



Davis et al.1 119851). and the spherical overdensity finder 
Lacev fc Cold Il994l ). Note that different algorithms may 



lead to differences of the o rder of 10-20% in th e mass func- 
tions (see the discussion in ljenkins et alj|200ll ). 

The UMF is related with F(m, z), the mass fraction in 
collapsed objects with masses greater or equal than m, as 



n(m, z) 



p dF(m, z) 
m dm ' 



(2) 



where p is the mean matter density of the Universe. It is also 
common to write the mass dependence in terms of a(m,z), 
where a 2 (m,z) denote the mass variance of the linear den- 
sity field extrapolated to the redshift z at which the haloes 
are defined, and on mass scale m oc ~pr 3 (the precise form of 
the m — r relation depends on the window function used). 
This variance can be obtained as 



a 2 (m,z) = ^l I \5 k \ 2 W 2 (kr)k 2 dk 



(3) 



where \8k\ 2 is the linear power spectrum of density fluc- 
tuations, and b(z) is the growth facto r of linear pertur- 
bations normalised to unity at z = (|Lahav et al.l Il99ll : 
ICarroll et aT]|l992h . The window function is usually taken 
to be a real-space top-hat filter, i.e. m = 47rpr 3 /3, and 
3 

W(x) = — (sinx — xcosx) (4) 

Using lner -1 as the mass variable, the UMF is also de- 
fined in some works as 

m dN(m, z) 



p diner -1 



(5) 
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where now / represents the mass fraction contained in col- 
lapsed objects per unit In a -1 , and N(m,z) is the number 
density of collapsed objects with masses above m, i.e. 



N(m y z) = / n{m' , z)dm! . (6) 

J m 

The / function is connected with the standard F(m, z) as 



F{a,z) 



f(<j', z)dln a' 



(7) 



The advantage of using eq. [5] is that the majority of 
the analytic (or semi-analytic) models for the UMF pre- 
dict a functional form for f(a, z) with no explicit depen- 
dence on redshift. In those cases, all the dependence in 
mass and redshift can be absorbed into a single variable, 
v = (S c /cr(m, z)) 2 , where S c is the threshold parameter usu- 
ally taken to be the extrapolated linear overdensity S\ of a 
spherical perturbation at the time it collapses. For the case 
of an Einstei n-de Sitter cosmology, we have 8 C = 1.686 (e.g 
|Peebleslll980l ). with a weak depende nce on the cosmological 
parameters (e.g lNavarro et aT]|l997T l . 

In this paper, we shall consider four of the different pre- 
scriptions for the UMF which were mentioned in the intro- 
duction, namely PS, ST, WA and BM. However, the method 
described here can be applied in principle to any other pre- 
scriptiorQ. The expressions for these UMF are the following: 

(i) The Press - Schechter mass- function 

i Press fc Schechterl 1 1974 iBond et all Il99ll ; lLacev fc Coiel 
1 19931 ) is given by 



/psO) 



2 5 C 



exp(- 



2a 2 ' 



(8) 



(ii) The empirical best-fit mass function of ST is given by 



/ ST (a) = A 



1 + 



a 

2S 2 



•exp( 



aS 2 , 
"2a* ■ 



(9) 



where A = 0.3222, a = 0. 7 07 an d p = 0.3. 

(iii) The IWarren et ail (|2006l ) UMF uses a modified ver- 
sion of the functional form proposed by ST which is fitted 
to the simulations, giving 



/wa(ct) = A(a a + b) exp( -) 



(10) 



with A = 0.7234, a = 1.625, b = 0.2538 and c = 1.1982. This 
expression provides a very good fit, over a mass range of 5 
orders of m agnitude, to the UMF obtained from simul ations. 

(iv) The iBetancort-Riio fc Montero-Dortal (|2006al 'l mass 
function does not depends only on v, but there is an addi- 
tional dependence on the mass which enters through certain 
local spectral index, which is characterised by a parameter 
represented by c(m). In this case, the UMF is given by 



/bm(c, c(m)) = 



<1Fbm (m) 
din a 



(11) 



where the mass fraction .Fbm(wi) is defined as 



1 As we will see in Sec. [5] in order to use the proposed normalisa- 
tion condition, it is required that the mass function should have 
appropriate asymptotic behaviours at high and low mas ses. For 
this reason, we can not use here the I Jenkins et aL | toll) UMF. 



F B M(m) = 
V(m) = 3 



■FsT,a=i(m) 
V(m) 



erfc 



exp(- 



exp(— c(m)u 2 ) 



1/2- 



(12) 



u du (13) 



and FsT,&=ii m ) is the ST mass fraction which corresponds 
to equation [9] but using a — 1. The expression for c(m) will 
be discussed in Section [4.31 f equations l24l 1271 and I28 |l . 

Given that the BM mass function is defined in terms of 
the mass fraction F(m), we find it useful to have analyti- 
cal expressions for the mass fractions of the different UMF. 
These formulae are presented in Appendix lAl 



3 THE CONDITIONAL MASS FUNCTION 

As discussed above, one can find several analytic expressions 
in the literature for the conditional mass function (CMF) of 
dark matter haloes, n c (m). The most widely used frame- 
work is the so called extended Press-Schechter (EPS), or 
excursio n set formalism, either in the contex t of spherical 
collapse <|Bond et al.lll99ll ; lLacev fc Colelll993l ) or in a more 
elaborated form, in the case of a constant barrier. According 
to the later formalism, the CMF can be estimated by con- 
sidering successive crossings of barriers of different heights. 
In the case of constant barriers this problem has an analytic 
solution, and it is possible to provide analytic expressions 
for both the UMF and the CMF. In particular, the CMF es- 
sentially has the same form as the UMF, but in the rescaled 
variables 



S c — So 
a -ctq, 



(14) 

where So and o"q are the linear density and the amplitude 
of the condition. Hereafter, we will refer to this equation as 
the standard re-scaling. Note that this is precisely the scaling 
which is obtained for the conditional probability distribution 
of the average density field within scale a 2 at a randomly 
chosen point within a region on scale a 2 with inn er linear 



density contrast So (see appendix A of lBower1ll991 ). Within 
this EPS formalism, the CMF would be given bjE 



n c ,EPs{m\So,<Jo) 



exp ■ 




(15) 



However, as discussed in the last section, the EPS is not the 
most accurate approximation to the halo abundance, being 
one of the reasons that it does only contains the physics of 

the sphe rical collapse. 

The ISheth fc Tormenl (|200l hereafter ST02) approxi- 
mation for the CMF further extends this EPS formalism by 
including the physics of the ellipsoidal collapse. This is done 

2 Note that according to this definition, the conditional mass 
function n c is defined as "lagrangian" , i.e. we explicitly use the 
mean density, p, when defining the volume element, and not the 
actual local density in the conditional region. Throughout this 
paper, we will always refer to lagrangian CMFs unless otherwise 
stated. 
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by us ing the moving barrier shape derived by ISheth et al.l 
l|20011 ). which is given by 



-Bec(c 



= y/a5c 



aS c (z) 



(16) 



with a = 0.707, (3 = 0.485 and a = 0.615. Note that with 
our notation, the redshift dependence on the growth fac- 
tor is contained inside a(m,z), so the residual dependence 
8 C = 8 c (z) is due to the considered cosmology. Unfortu- 
nately, for this problem there is no analytic formulae for 
the first-crossing distribution. Nevertheless, ST02 proposed 
a simple analytic expression for the CMF, based on the di- 
rect replacement of the barrier shape in their generic ap- 
proximate expression for the first-crossing distribution. This 
proposal reasonably reproduces the results from numerical 
simulations for the mass function of the progenitors of haloes 
in a given mass range todaj0. The ST02 proposal can not be 
directly applied to obtain the CMF at a fixed redshift z and 
for the conditioning we are considering in this paper. How- 
ever, we use here a natural extension of the ST02 proposal, 
which consists in the replacement B(a 2 , z) — » B(a 2 , z) — (Jo- 
in that case, the explicit expression for this CMF is 



:,ST02(m) = - — 7— 

V 7r / m (a 



\T(a 2 \a 2 )\a 



x exp ■ 



(a 2 - a 2 )3/2 

[B cc {a 2 ,z)-5 ] 2 
2(a 2 - a 2 ) 



da 
dm 



where 



T(aVo) 



\ - {<J 2 -cj 2 ) n d n lB ec (a 2 ,z)~5 ] 



n=0 



d{a 2 Y 



(17) 



(18) 



As we show below, this prescription provides better results 
than the EPS for underdense regions (So < 0), although 
still shows a discrepancy with the data from simulations, 
underestimating the number of haloes at low masses. 



4 EXTENDING THE UMF 

Here we present our method to build the CMF of dark mat- 
ter haloes. An important formal difference with respect to 
the expressions discussed in the last section is the explicit in- 
troduction of the dependence in the q variable, the distance 
from the center of the condition to the point at which we 
evaluate the CMF. Thus, in order to fully describe the CMF, 
we have to specify three parameters, namely Q (the radius 
of the condition), So (the linear density within Q) and q. 
Thus, 



4.1 The method 

The method pr oposed here is an extended version of the 
one presented in lPatiri et alj l|2006l . hereafter PBP06), which 
was based on a set of assumptions similar to all other pre- 
viously described CMFs in the literature. For clarity, we 
follow the same notation as in PBP06. The basic idea of the 
method is to build the CMF as an analytic extension of any 
of the existing versions of the UMF. To this end, we make 
two basic assumptions: 

(i) For those scales m which are much smaller than the 
scale of the condition m(Q), the constrained field should be- 
have "locally" as an isotropic uniform Gaussian field with 
a re-scaled mean and variance (or power spectra). This as- 
sumption is implicitly done in all the existing derivations of 
the UMF (although for the average value of the field). For 
a discussion on its validity, see PBP06. 

(ii) We assume that the UMF can be derived from a 
"rigid" barrier, B(m), in the sense that if we have a crossing 
of this barrier, then we will have collapse with probability 1 
(note that for this argument B(m) may have some depen- 
dence in the mass) . In the general case of ellipsoidal collapse 
the barrier should have certain "width" , connected with the 
fact that the evolution of an ellipsoidal perturbation is de- 
termined by three parameters (the three eigenvalues of the 
deformation tensor) . Thus, for a given mass (or equivalently, 
a given a) there is no a rigid height of the barrier, but we 
have a probability distribution of finite width around a cer - 
tain value of S\. However, as discussed in ISheth et all (|200ll ). 
using an average (rigid) shape for the barrier for treating the 
ellipsoidal collapse gives good results for the UMF. 

After these two assumptions, and neglecting any depen- 
dence of the shape of the barrier which is not captured by 
the ratio S c /a, it follows that obtaining the CMF just im- 
plies a vertical displacement of the barrier (without changing 
its shape). This displacement is obtained as a re-scaling in 
which the variables S c and a(m, z) are substituted by their 
corresponding values in the local Gaussian field. Note that 
a more detailed treatment of the ellipticity of the collapse 
would require to compute the exact shape of the barrier after 
each re-scaling, implying an integration over the distribution 
function for the three eigenvalues of the deformation tensor. 

This basic idea of re-scaling the UMF is identical to the 
one adopted in PBP06. However, the overall normalisation of 
the resulting CMF was not treated in detail in that work, as 
we explain below. Here, we extend the method by obtaining 
this normalisation as it is described in Sec.JSJ and we present 
the full explanation of the approximate treatment of the 
normalisation used in PBP06. 



n c = n c (rn, z\Q, So, q). 

Note that the quantities Q and q denote the Lagrangian 
radius; for Eulerian radius we will use R and r, respectively. 



3 E.g. if we consider the conditional mass function of objects 
with variance <ti at redshift z\, given that we have a variance 
(To at an earlier redshift 20, their proposal for the re-scaling, to 
be inserted in their equation for the first-crossing distribution, is 
B(a 2 ,z) -> B{a\,zi) - B(a 2 ,z ). 



4.2 Formalism 

The relevant statistical quantity is the conditional probabil- 
ity distribution, P(S2\S\ , q, Q) for the linear density pertur- 
bation ($2, on scale Q2, at a distance q from the center of a 
sphere of radius Q (the condition) with mean inner linear 
density fluctuation Si. For a Gaussian field, this is 



P(S 2 \Si,q,Q) = (2n) 



-1/2/ 



02 



4y 



1/2 
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exp i 



(19) 



where a\ =< 5j >= a 2 (Q), a\ =< of >= cr 2 (Q 2 ), and 
(Ti2 = 0-12(9, Q,Qi) = 

1/ \ /> + oo 

lZJ ' |<5 fe | 2 W / (fcQ)VF(fcg2)sin(fcg)fcdfc. 



(20) 



If one compares the one-point statistic for 5 2 with the 
unconstrained case: 



P(5 2 ) = G(5 2 ;a 2 ) = (2 7 r)- 1/2 a 2 - 1 exp ( - ±i|) 



(21) 



we see that the conditional case behaves in the same way, 
but with re-scaled variables S[ and a' given by 



51 



61 



<7\1 



5 2 -6iD(q,Q,Q 2 ) 



(22) 



The BM mass function shows an additional depen- 
dence on mass apart from the standard v dependence, 
which is included in the c(m) coefficie nt, defined as (see 
iBetancort-Riio fc Montero-Dort j|2006bl ) 



c(m(Q)) = -ln[D(Q,Q,Q)] 



(27) 



In order to re-scale this UMF, we use the fact that c(m) ~ 
^d\na/d\nr, so the corresponding scaling for the c(m) term 
will be 



, ldlna' 
c(m) — > c (m) « — - 



1 do'- 



= c(m)l- 



(28) 



where a' is defined in equation 1251 

Finally, we note that by definition, all the UMFs verify 
that f(y) — if v < 0. This has to be explicitly taken into 
account because D(q,Q,Q 2 ) can be larger than 1 for some 
values of q (for example, at q = 0), and thus, one may find 
v' = 5'c /a' 2 < even if 6 C > 61. 



and 



o-'(m) = (crl - ^) V2 = (a 2 2 - D(q, Q, Q 2 ) 2 o?) 1 2 5 ! 
where we have introduced the functiorQ 

0-12(9, Q,Qi) 



1/2 



D(q,Q,Q 2 ) = 



(24) 



In Appendix [B] we present some useful fits to this function 
in several asymptotic cases. One of the important properties 
that can be used to speed up the computations with this 
function is that if Q 2 -C Q, then D(q,Q,Q 2 ) is practically 
independent of Q 2 . 



4.3 Re-scaling the UMF 

From the previous discussion, it is natural to adopt the fol- 
lowing change of variables in order to re-scale locally a given 
UMF to obtain the corresponding CMF: 



8' c = 5 c -D(q,Q,Q2)5i 
(a') 2 =a 2 -D(q,Q,Q 2 ) 2 o 2 



(25) 



In principle, we can use any of the aforementioned expres- 
sions for the UMF to build the CMF. However, there are 
some details which are important to discuss. 

In the first two cases considered for the UMF (PS and 
ST), the UMF can be explicitly written in terms of the vari- 
able v = (S c /a(m, z)) 2 . Thus, re-scaling these mass func- 
tions just implies the change of variables given in eq. 1251 
However, for the WA mass function only the dependence on 
a 2 is explicit shown, so to use this UMF (or any other nu- 
merical fit) we proceed as follows. We shall assume that all 
the dependence on the mass can be absorbed in the variable 
v — (1.686/ct) 2 , and we will re-scale this variable accord- 
ing to eq. 1251 In practise, this implies making the following 
change in eq. [10] 



a 2 - D(q,Q,Q 2 ) cr 1 



1.686 - D(q, Q, Q 2 )5i 



(26) 



4.4 Examples and comparison with other 
prescriptions 

Fig. [T] shows two particular examples of the CMF, built 
in this case from the ST mass function, and using Q = 
8 h _1 Mpc (equivalent to w 1.8 x 10 14 IT 1 M©). We consider 
the case of an underdense region (Si = —4), and an over- 
dense one (Si = 0.5). As discussed above, our prescription 
provides the CMF as a function of the position within the 
condition, so the figure illustrates the range of variation of 
n c with the radial coordinate q. For the case of (5i = —4, the 
radial dependence is strong for those masses which are a sig- 
nificant fraction of the total mass of the conditioning region, 
being the number density significantly smaller at the center 
of the condition. For the case of Si — 1, the radial depen- 
dence is important practically at all masses. In particular, 
at low masses the number density is significantly smaller at 
the center of the condition than at the boundary. One can 
easily understand this result in this way: near the center, 
D(q) w 1.3, which makes the value of the local density very 
close to the threshold. In that case, most of the mass will 
be collapsed in large objects (with masses of the order of a 
significant fraction of m*, with o(nu) = 1), and low mass 
object will be less abundant in number when compared with 
the border. 

To compare our prescription with results existing in the 
literature, we need to obtain the "volume average" CMF 
over the condition, nf s . This function is easily derived as 



n^(m,z\Q,Si) 



n c (m, z\Q, 81, q)q 2 dq 



(29) 



Note that this D(q) function is independent on redshift. 



For such comparison, we also show in Fig. [T] this average 
function, which, as a consequence of the radial dependence 
commented above, will be much more close in shape and 
amplitude to the CMF with values of q in the vicinity of Q. 

Fig. [2] compares these two average CMFs with other 
existing prescriptions. In particular, we first want to an- 
swer the following question: does our averaged CMF coin- 
cides with the CMF that would be obtained using the "stan- 
dard rescaling"? In this context, we consider as "standard 
re-scaling" the one presented in equation (|14p . but applied to 
any of the aforementioned UMFs and not only the PS UMF. 
The interesting result shown in Fig. [2] is that, for a given 
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Q= 8.0 h"' Mpc, 6, = -4.0 Q= 8.0 h" 1 Mpc, 6, = 1.0 




10'° 10" 10 12 10' 3 10 10 10" 10 12 10' 3 

m [rf 1 MJ m [V 1 MJ 



Figure 1. Radial variation of the CMF. Two examples of the CMF, n c (m, z\Q,5\,q), are shown for the case of an underdense region (<5i = 
—4, left panel) and an overdense region (<5i = 1.0, right panel), for a condition of Q = 8 h^ 1 Mpc (equivalent to Rj 1.8 X 10 14 h" 1 Mq). 
The CMF is obtained using the re-scaling presented in equation 1251 and the reference UMF is taken to be the ST mass function. Both 
panels present the two extreme cases of q = (solid line) and q = Q (dashed line), showing the range of variation of the CMF within the 
condition. For comparison, we also present (dotted line) the average CMF over the condition (which is computed using equation [29} . 




Figure 2. Comparison of the average CMF from four different prescriptions. The first two cases correspond to the standard rescaling 
(see text for details) , and to the average CMF (n^ VE ) computed as described in equation 1291 Both of them are computed using the 
ST UMF. The other two prescriptions correspond to the EPS (eg. I15t and to the ST02 (eg. 1171 . As in Fig. [T] the left panel refers to 
an underdense region (5i = —4), while the right panel corresponds to an overdense one (<5i = 1.0), both of them for a condition of 
Q = 8 h" 1 Mpc. 



UMF (e.g. the ST in this case), the corresponding average 
CMF (computed with our formalism using equation I29p , co- 
incides with great accuracy with that CMF obtained with 
the "standard re-scaling" , but only for those masses much 
smaller than the condition. However, and specially for the 
case of underdense regions, when we consider larger values 
of the mass (in this particular value of the linear density, 
for masses > 1CP 2 m,), then the "standard" computation 
differs from the exact result. 



different results. For low masses, the dependence of the mass 
function on the local linear density fluctuation (whose mean 
value at q is 5\D(q, Q, Q2)) is close to linear, and thus the 
processes of rescaling and averaging almost commute. How- 
ever, for higher masses, non-linear terms are important and 
this two processes do not commute. 



One can understand this result in the following way. 
Within our approach, we are rescaling the mass function lo- 
cally, according to the local Gaussian field, and after that, 
we are averaging over the whole volume. However, in the 
"standard case" , the volume- averaged re-scaling of the den- 
sity field is used as an ansatz for the prescription to rescale 
the UMF. In principle, these two procedures could produce 



For comparison, Fig. [2] also presents the CMF obtained 
from EPS (eq. [TBJ and ST02 (eq. (TTJl prescriptions. All the 
different mass functions give similar results for the overdense 
case, while showing significant discrepancies for the under- 
dense case. 
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5 A NORMALISATION CONDITION FOR 
THE CMF 

In order to assess the quality (and the validity) of a given 
expression for the CMF, we propose here to use a "normali- 
sation condition" that connects the CMF with the uncondi- 
tional one. This condition must be satisfied by any CMF in 
and reads 



Q = 8 h"' Mpc 



n{m) = / n c (m\Si,Q,q — 0)G(Si;ai)dSi, m < m(Q)(30) 

J — oc 

here G(<5i;<7i) stands for a normalised gaussian with mean 
Si and variance ai, and we have omitted the redshift depen- 
dence in both sides for simplicity. Note that, by construc- 
tion, there can not be masses larger than m(Q) within the 
conditioning region. It is also important to stress that in 
this equation, the lagrangian CMF is the one that has to be 
used. 

Equation[30]quantifies the fact that when one integrates 
the CMF for all possible values of the linear density within 
the condition (Si), then one should recover the UMF. By 
definition of the CMF, it is clear that this equation is exact 
in the limit eri — > (in that case, G(Si;ai) — > S^(5i), 
where <5 ( - D - ) is the Dirac delta- function). 

For finite values of oi(Q), it should also be satisfied 
to a good approximation for large values of Q (i.e. in the 
linear regime), and for those masses m much smaller than 
m(Q). The reason for this is related with the way we have 
set the upper limit of integration. In principle, one would 
recover the true UMF if the integration is carried out for all 
possible values of Si. However, the CMF will not be defined 
for values of Si above the barrier, because in that case the 
whole conditioning region would have collapsed, and only 
masses larger than m(Q) may exist (i.e. it would be part of 
a larger object). Therefore, under the assumption of a rigid 
barrier, the upper limit of the integral in eg. 1301 should be 
given by this barrier, i.e. 8 C in the case of spherical collapse, 
or by B cc (ai,z) (which wa s defined in equat ion \T6$ for the 
ellipsoidal collapse barrier (|Sheth et al.ll200ll ). 

As discussed above, this abrupt way of truncating the 
integral is based on the assumption of rigid barrier, and it 
might fail for values of Si close to the boundary (<5 C ). For 
this reason, and in order to normalise the CMF, we will 
focus on large values of Q for the computations (<7i < 1). 
In this way, the gaussian function in the rhs of eq. [301 is 
narrow, and the details of how the truncation is done are not 
important, while still eq. [30] imposes a significant constraint 
on the shape of n c . 

Fig. [3]presents the relative difference between both sides 
of Eq. 1301 defined as the ratio R(m) of the right- to the left- 
hand sides of that equation, i.e. 

f X n c (m\Si,Q,q = 0)G(Si; ai)dSi 

R(m) = ^22 (31) 

n(m) 

where X is equal to B ec (ai, z) for those cases in which the 
UMF includes the physics of the ellipsoidal collapse, and 
X — S c for the PS case. For illustration, we consider three 
cases for the mass function to be used with the standard 
re-scaling procedure, namely ST, WA and BM, and the con- 
dition R c = 8h _1 Mpc (or 1.79 x 10 14 h" 1 M Q ). The main 
conclusion from that figure is that none of the considered 
mass functions satisfy the normalisation condition, showing 




Figure 3. Relative error in the normalisation equation 1301 
for the CMF. We consider the case of R c = 8 Mpc/h. This 
percentage error is shown as a f unction of the mass m for 
three different mass functions: ST dSheth fc Tormer]|l999l). BM 
[|Betancort-Riio fc Montero-Dortall2006ah and WA l lWarren et al.l 
l2006t ). It can be seen that for the ST and the WA cases, the rela- 
tive error is larger than 10% for high masses. For comparison, it 
is also shown the relative error for the case of the EPS and ST02 
prescriptions. 



deviations in some cases larger than 20%, specially for high 
masses. If we repeat this computation for different values 
of the condition, we find that for a given mass, the relative 
error becomes larger for smaller conditions, or equivalently, 
for larger values of <ti (for example, for R c = 6h _1 Mpc 
and the ST mass function, the relative error becomes al- 
ready of the order of 5% for m « 10 11 h^ 1 M©). This fact 
shows that the re-scaling proposed in that case to build the 
CMF can not be exact. For comparison, we also present the 
corresponding R(m) function for the EPS and ST02 pre- 
scriptions, which can be computed because the normaliza- 
tion condition (eq. 130 also applies to the averaged CMF. 
We note that the EPS is practically normalised, although it 
does not includes the physics of the ellipsoidal collapse. On 
the other hand, the ST02 satisfies the normalization condi- 
tion with high accuracy (better than approx. 5 per cent), 
although as it is shown below in section [7] it does not repro- 
duces the numerical results in simulated voids. 



6 EXTENSION OF THE METHOD 

The incorrect normalisation of the CMFs built using the 
standard re-scaling as described in the last section imply 
that at least one of the assumptions made in the prescrip- 
tion for re-scaling the UMF is incorrect. Given that the hy- 
pothesis of a rigid barrier has provided very good results for 
the case of UMFs, and that the hypothesis of local isotropy 
should be satisfied to great accuracy, the only remaining 
assumption seems to be the adopted expression for the re- 
scaling. In PBP06, the question of the normalisation was 
handled in a simple manner: the re-scaled mass function was 
simply divided by R(m). This "first-order approximation" 
procedure provided good results, but this a posteriori renor- 
malisation can not be entirely correct. In principle, if one has 
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Table 1. Values for the 5 c g parameter that satisfy the normal- 
isation condition (Eq. I30|l for the different UMFs considered in 
this paper. These values have been obtained from a condition of 
6 h _1 Mpc. For each one of the UMF presented in the first col- 
umn, the second column shows the interval in which the error 
in the normalisation condition is smaller than 5%, at the largest 
considered mass (Q V oid/30). Third column shows the "optimum" 
5 e g which we finally adopt, obtained as the value that gives the 
smallest percentage error in the normalisation condition at the 
highest considered mass (Q V oid/30). 



UMF 


<5 c ff 


Sett 




[Err < 5%] 


(Adopted) 


PS 


[1.31,1.53] 


1.42 


ST 


[1.20,1.42] 


1.25 


WA 


[1.02,1.19] 


1.10 


BM 


[1.37,1.64] 


1.50 



a correct re-scaling procedure, the normalisation condition 
must follow automatically. 

A complete study would require a detailed derivation 
of the exact shape of the ellipsoidal barrier in a conditioned 
environment, which is beyond the scope of this work. In- 
stead of that, we propose here a simple modification of the 
re-scaling law in order to fulfil the normalisation condition. 
Our choice to propose a modification of eq. [5S]is to assume 
that the linear density which has to be re-scaled is not S c , 
but a certain parameter 8 e a- We stress that this is simply 
an ansatz to absorb possible errors implicit in the re-scaling 
procedure given in equations [14] and 1251 so as to render a 
modified prescription which would satisfy the normalisation 
condition with much higher accuracy. This means that we 
should not attach much meaning to the values of <5 e ff ob- 
tained. However, as we note latter, the fact that a particu- 
lar UMF satisfies the normalisation condition with a value 
of <5 G fi of the order of 1.6, may be considered as a "natural" 
result and a strong point in favor of that UMF. 

If the variable that has to be rescaled is S e s instead of 
<5 C , then the new recipe would be 

8 c ^S' c = 6 c y--D(q,Q,Q 2 )p^j (32) 

or equivalently, the total re-scaling for the variance will be 

2 _ a 2 ~D(q,Q,Q 2 ) 2 ai 
a ~* (l-D(q,Q,Q 2 )5 1 /S cS r [ ' 



6.1 Calibrating the 8 c g parameter 

Using this new prescription (equations 1321 and I33 [) to build 
the CMF, we can solve the normalisation condition (equa- 
tion [301) to find a solution for 8 e g, which in principle would 
be a function of mass. Given that the different expressions 
of the UMF are not exact, but they have numerical uncer- 
tainties (typically of the order of few percent), we decided 
to find, for each mass, the range of 8 e g values that fulfil the 
normalisation condition within a certain accuracy. Then, we 
look for a single value of S e g which satisfies eq. [30] with high 
accuracy for all masses. 

Figure[4]presents the allowed region for 8 e g as a function 



of mass, which fulfil the normalisation condition for a cer- 
tain UMF with an accuracy better than 10, 5 and 3 per cent, 
respectively. For this figure, the case of Q — 6 h" 1 Mpc was 
considered, and the upper limit of the integral in equation l30l 
was taken to be the barrier, B ec (a 2 , z). As one would expect, 
for a given UMF the strongest constraint on 8 e g is obtained 
for high masses, because at those masses the differences be- 
tween the CMF and the UMF become more important. If 
our ansatz were the correct scaling, then one should be able 
to find a value of S g for which the error in the normalisation 
would be zero for all masses. This is not the case, for exam- 
ple, for the PS UMF, where the "allowed region" changes 
with mass. The reason for this could be that this is the only 
UMF which does not take into account the physics of the 
ellipsoidal collapse. However, the important point is that in 
all the other cases (ST, BM and WA) it is possible to select 
a value for 8 e g which makes the error in the normalisation 
condition to be smaller than ~ 5 — 10% for all masses. 

Another important issue is the dependence of the <5 e ff 
with the size of the condition, Q. If our ansatz were the cor- 
rect scaling, then this parameter should not depend on Q. 
We have checked that the 5 e ff values do not show a strong 
dependence on the chosen size for the condition, in the sense 
that the range of 8 e g which fulfils the normalisation equa- 
tion with accuracy better than 5% is always overlapping 
when changing the size of the condition between 5 h _1 Mpc 
and 15 h~ Mpc. Indeed, choosing a value too high for the 
condition does not provide a good constraint on <5 e ff, because 
in that limit the conditioning fluctuation is very small and 
it is not imposing strong constraints on the exact shape of 
the CMF. We illustrate this fact in Figure[5] where we show 
the allowed regions for the case of the ST UMF, but consid- 
ering now the condition Q = 10 h~ Mpc. For this reason, 
we decided to adopt as reference values for this paper those 
obtained with the case Q = 6 h _1 Mpc. 

The adopted values of <5 e ff for the computations in the 
rest of this paper are presented in Table [1] Our criteria was 
to take as "adopted" 5 e g value the one giving the lowest 
relative error at the highest considered mass in each case, 
which in our computations was taken to be m(Q)/30. How- 
ever, we note that as regards to the normalisation condition, 
any of the values quoted in the 5 per cent range could be 
adopted. 

6.2 Examples of the CMF with the 8 e g parameter 
included 

The values of S e g obtained for different fits to the UMF 
are somewhat different, but for each particular case, they 
lead to similar shapes for the CMF. We illustrate this in 
figure [5] where we present the average CMF for the two 
cases considered before (Q — 8 h" 1 Mpc, with Si = —4 or 
Si — 1), and using several UMFs with their corresponding 
"optimum" 8 e g values. Given that the normalisation condi- 
tion is fulfilled by all the different CMFs with a precision 
of ~ 5%, one would expect differences of this order between 
them. However, this is not the case. For the low-density case 
(5i = —4), there are significant differences at high masses, 
although the agreement at low masses is very good. And for 
the high density one (Si = 1), the discrepancies are larger 
than that even at low masses. 

To conclude this subsection, we revisit again the discus- 
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PS mass function (Q = 6 h 1 Mpc) 



ST mass function (Q = 6 h 1 Mpc) 




WA mass function (Q = 6 h 1 Mpc ) 

i i n i -n i i i i ^ 

10% 
5% 
3% 





m [h" M»] 

BM mass function (Q = 6 hT 1 Mpc) 




Figure 4. Ranges of the <5 e ff parameter which make the normalisation condition to be fulfilled with certain precision (10, 5 and 3 per 
cent, respectively), as a function of the mass. We consider four different prescriptions for the mass function: PS (upper left), ST (upper 
right), WA (lower left) and BM (lower right). All computations have been done for the case of Q = 6 h — 1 Mpc. 
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Figure 6. Comparison of the average CMF (nj vg ) obtained within our formalism (i.e. using the scaling given in eqs, I32I and |33> . but 
using different UMFs. We consider the ST, BM and WA UMFs, with their corresponding <5 e ff value taken from Table[T] and for the same 
particular cases discussed before in the paper (Q = 8 h" 1 Mpc, and either 5i = -4 or Si = 1.0). 
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Figure 7. Comparison of the average CMF from four different prescriptions. This figure is equivalent to Fig. [2] but now we use the 
optimum 5 ff instead of 8 C when plotting the CMF. The first two prescriptions correspond to the standard re-scaling (eg. I14|l . and to the 
method proposed in this work (eq.[32j, both using the ST UMF. so 8 cB = 1.25 according to Table[T] The EPS and the ST02 CMFs are 
also shown for comparison. As in Fig. [2] we consider a condition of Q = 8 h Mpc, and &\ = —4 (left panel) or <5i = 1.0 (right panel). 



ST mass function (Q = 10 h ' Mpc) 




10" 10 ,z 
m [h~' Mj 



Figure 5. Same as figure|4] but for the case of Q = 10 h — 1 Mpc. 
Although the allowed region at the level of 5% is compatible with 
the one obtained for Q = 8 h _1 Mpc in Fig. [4] in this case it 
is much larger, and thus less restrictive if one is interested in 
constraining the value of 5 e jf. 



sion about the difference between our averaged CMF, and 
the CMF obtained with the standard (average) re-scaling of 
the UMF, but now using 5 c g instead of 5 C as the variable to 
re-scale. We have already shown in Sect. l4.4l that for the case 
of 5 e ff = 5 C , the standard re-scaling provides good results for 
low masses (see Fig. [5J). As a example, we now consider in 
Figure [7]the case of 4ft = 1-25 for the ST UMF. The result 
for the under-dense region (Si = —4) is qualitatively similar 
to the one obtained before for the case of <5 e ff = <5 C - Indeed, 
the discrepancies between the averaged CMF and the CMF 
with the average standard re-scaling appear approximately 
at the same mass scale (m ~ 10 -2 m*). Thus, the "average" 
re-scaling still provides good results for those masses much 
smaller than the size of the conditioning region. However, 
this is not the case for the over-dense (<5i = 1) region, in 
which even at low masses, there is a discrepancy with the ex- 
act result. The reason for that is that in this case, the region 



near the center has surpassed the critical value of <5 e ff (note 
that 5eff = F25 for the ST UMF, and that D(q = 0) w 1.3), 
so in practice that region is not contributing to the average. 
As a consequence, the CMF built with standard re-scaling 
slightly overestimates the number of low-mass objects. 

Finally, it is interesting to compare these results with 
the EPS and ST02 prescriptions, which are also presented 
in Fig. [7] For the under-dense region (<5i = —4), our final 
version of the CMF is now in closer agreement with the 
ST02, although there is still a systematic offset at low masses 
between the two prescriptions of the order of 30-40%. The 
discrepancy with the EPS is much larger. As can be seen 
from the numerical results presented in the next section, this 
later prescription significantly over-estimates the number of 
haloes at low masses. For the over-dense region (Si = 1), 
both EPS and ST02 are systematically above our final result. 



7 THE CONDITIONAL MASS FUNCTION IN 
VOIDS 

One of our main motivations for the development of an ac- 
curate prediction for the CMF is the theoretical study of 
the statistic of voids. PBP06 developed a general analytical 
procedure for computing the number density of voids with 
radius above a given value, which was applied successfully to 
the description of the statistics of voids found in numerical 
simulations. 

Their formalism is based on a detailed study of the num- 
ber density of non-overlapping empty spheres with radius r. 
PBP06 provided an analytical expression relating this quan- 
tity with the void probability function (VPF), Po(r) (i.e. 
the probability that a randomly placed sphere of radius r is 
empty). Po(r) can be obtained from two basic ingredients, 
namely the probability distribution function for the values of 
the density contrast, S, within a randomly chosen sphere of 
radius r; and 8n(Si), the mean fractional fluctuation within 
r of the number density of the objects defining the void as 
a function of the linear fractional density fluctuation within 
r. In detail, the final expression is 
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Po(r) 



-rTV[l+« N («i)] 



P(Si\r)dS l 



(34) 



In this section we will focus in the 8n(8i) function, be- 
cause it encodes all the information about the CMF. This 
function can be written as 



l + S N (m,Q,Si) = [l + S ns (m,Q,Si)][l + S(6i)] 



(35) 



where 8 ns accounts for the "statistical fluctuation" (i.e. the 
clustering of the protohaloes in the initial conditions be- 
fore they move with mass), and S(Si) is the actual fractional 
mass density fluctuation as a function of the linear value. As 
shown in PBP06, the S ns term can be related to the CMF 
in the following way 



1 + 8 nB (m, Q,S\) = 



N(m) 



-p^ I N c (m\Q,S 1 ,q)q dq 




(36) 



where we have introduced the number density of collapsed 
objects with masses above m for the conditional case, which 
is given by 



N c (m\Q,8 h q) 



n c (m'\Q, 5i, q)dm! 



(37) 



We can use now our proposed formalism to produce 
accurate computations for the S n s(m, Q, Si) function which 
can be used when applying the aforementioned formalism. 
As proposed in PBP06, equation [36] can be well-described 
with a function of the shape 



1 + 8 nB (m, Q, Si) = A(m, Q)e 



-i>(m,Q)«f 



(38) 



We have checked that this fitting formulae provides a reason- 
able fit for those masses m > 3 x 1CF 3 nu (with cr(nu) = 1). 
Detailed computations of the A(m, Q) and b(m, Q) functions 
for the cosmology considered in this paper (as = 0.9 and 
T = 0.21) are given in appendix [C] The dependence on cos- 
mology of the different coefficients of the fits is presented in 
appendix iDl 

Finally, as an illustration of our method, we compare 
the CMFs of haloes within voids with the results o f numeri- 
cal simulations presented in lGottlober et ail (|2003l ). Figure 5 
in that paper presents the mass function of haloes in five sim- 
ulated voids for a cosmology with fl m = 0.3. Two of them 
correspond to voids with radius i? vo id = 10h~ Mpc and 
mean density (in units of the critical density) of 0.03 (i.e. 
the density contrast is S = 0.03/0.3— 1 = —0.900); while the 
other three correspond to voids of radius -Rvoid = 8h _1 Mpc 
and mean density 0.04 (or equivalently, 8 — —0.867). In 
order to convert the mean density values ( 8) into linear den- 
sity o nes, we use the expression given by ISheth fc Tormenl 



1.68647 



1.68647 



1.35 1.12431 



(1+5)2/3 (1 + 5)1/2 



0.78785 

(1+5)0.58661 



(39) 



For the two considered voids, we obtain 8\ = —5.092 and 
5\ = —3.995, respectively. 

Figures [8l [9l and 1 1 01 compare the numerical results with 
our prediction for the CMF, using three different choices for 



the UMF, namely ST, BM and WA. In these figures, it is 
represented the accumulated Eulerian mass function aver- 
aged within the sphere, which in our formalism is computed 



N CJS (m\Q,8i) = (l + 5 n3 (m,Q,8i))(l + 8(6i))N(m) 



(40) 



where we use the S ns function computed as in equation 1361 
All three cases are in excellent agreement with the sim- 
ulations. Note that considering the case of 8 e g — S c m 1.69 
overestimates the mass function, specially in the high-mass 
tail. Note also that choosing a value for 8 e B which is 0.1 be- 
low the optimum value, already shows differences, specially 
in the first set of simulations (-Rvoid = 10 h" 1 Mpc). In ad- 
dition, the ST02 prescription systematically underestimates 
these numerical results by approximately 40%, specially at 
low masses where the numerical uncertainty is much smaller. 
The EPS prescription significantly over-estimates these re- 
sults (by practically a factor of two), so it is not shown 
in these figures. Finally, we stress that there are no free 
parameters in this approach. The value for 5 e g is uniquely 
determined by the normalisation condition. 



8 DISCUSSION AND CONCLUSIONS 

We have developed a procedure to obtain from any analyt- 
ical expression for the UMF (either theoretically motivated 
or numerical fit) a CMF which is in very good agreement 
with simulations both for underdense and overdense condi- 
tioning regions. 

To do this, we first considered the standard re-scaling 
(see eq. I14p in the theoretically motivated expressions for 
the UMF (those where 8 C appears explicitly), and pointed 
out that to improve the results the re-scaling must be im- 
plemented locally. We showed that for underdense regions, 
this leads to a non- negligible change for masses of the or- 
der of m*/100, while for overdense regions, the difference is 
significant only at higher masses. 

Next, we showed how to rescale (locally) any UMF even 
when <5 C does not appear explicitly ( equations 125126] ) . 

We then pointed out that those re-scalings are not fully 
consistent, because the UMF is not exactly recovered after 
integrating over all possible values of the condition. That is, 
the "normalisation condition" (eq. I30p is not satisfied. The 
discrepancy between both members of this equation may be 
above 20 % for the higher masses, and this is much larger 
than the precision at which we know the UMF from numer- 
ical simulations (~ 5 %). 

To solve this problem, we modified the rescaling proce- 
dure replacing <5 C by a quantity S c s to be obtained for ev- 
ery UMF by solving the normalisation condition, which now 
becomes an equation for 8 e g. Using these values of <5 e ff, all 
UMFs render very similar CMFs for both overdense and un- 
derdense regions, while for the former case there is a larger 
scatter in the mass functions. 

Based on these results, we also present (see ap- 
pendix|D}, an accurate fit for the (accumulated) conditional 
mass function, which depends on redshift z, and the cosmol- 
ogy (ag and V), which can be used in a variety of problems 
in large scale structure studies. 

Finally, we point out that based on the derived val- 
ues for <5 e ff, there is a preference for the BM UMF against 
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Figure 8. A comparison of the CMF of haloes in voids with the numerical results of lGottlob er et al, N2003I '). Left panel shows, in dotted 
lines, two simulated voids with -R vo id = 10 h — 1 Mpc (i.e. Q V oid = ^void(l + <S) 1//3 = 4.6 h _1 Mpc) and <5i = —5.092. Right panel shows 
three simulated voids with i? V oid = 8 h _1 Mpc (i.e. Q V oid = 4.1 h — 1 Mpc) and <5i = —3.995. In both panels, the other lines show our 
prescription for the CMF when implemented from the ST UMF. We consider three cases for <5 e fj , namely 1.25 (solid), 1.686 (dashed), 
and 1.15 (dot-dashed). For comparison, the ST02 mass function is also plotted. It is seen that using the standard value of 1.686 does 
not reproduce the numerical values. However, using the <5 e ff value which fulfils the normalisation condition, we correctly reproduce the 
amplitude and shape of the numerical results. 




Figure 9. Same as figure[8] but using the BM UMF as a reference. The results are similar to the previous case. Note that for this mass 
function, the value of S e g that reproduces the numerical results is closer to 1.686. 



the rest of the UMFs. Different fits lead to somewhat 
different values of <5 e ff, but all of them are incompati- 
ble with theoretically motivated values which must sat- 
isfy 5 c g > 1.6. However, if we apply locally the re- 
sc aling procedure that follows from the approach used 
bv lBetancort-Riio fc Montero-Dortal (|2006al ) , which involves 
re-scaling certain local spectral constant in addition to the 
mean and variance of the density field, the resulting CMF 
is also in very good agreement with the simulations but 
now the value of S e g is theoretically "acceptable" (note that 
the range of 5% is the only one which is compatible with 
5 e ff ~ 1.6). This suggests that this last rescaling is the "cor- 
rect one" , although most of the results in this work may be 
equally well obtained with any other prescription. 



ACKNOWLEDGMENTS 

We thank Stefan Gottlober for kindly providing us in elec- 
tronic form the statistics of voids from the simulations de- 
scribed in Gottlober et al. (2003). 



REFERENCES 

Betancort-Rijo, J. E., & Montero-Dorta, A. D. 2006a, ApJ 
Lett, 650, 95. 

Betancort-Rijo, J. E., & Montero-Dorta, A. D. 2006b, ApJ 

Lett, 653, 77. 
Bond, J. R., & Efstathiou, G. 1984, ApJ, 285, L45 
Bond, J. R., Cole, S., Efstathiou, G, & Kaiser, N. 1991, 

ApJ, 379, 440 
Bond, J. R., & Myers, S. T. 1996, ApJS, 103, 1 
Bower, R. G. 1991, MNRAS, 248, 332 



CMF of dark matter haloes 13 




Figure 10. Same as figures [8] and [9] but using the WA UMF as a reference. Note that for this mass function, the value of <5 e ff = 1.69 
clearly over-estimates the number of objects at all masses, specially in the left panel. 



Carroll, S. M., Press, W. H., & Turner, E. L. 1992, 

ARA&A, 30, 499 
Davis, M., Efstathiou, G., Frenk, C. S., & White, S. D. M. 

1985, ApJ, 292, 371 
Efstathiou, G., Frenk, C. S., White, S. D. M., & Davis, M. 

1988, MNRAS, 235, 715 
Gottlober, S., Lokas, E. L., Klypin, A., & Hoffman, Y. 2003, 

MNRAS, 344, 715 
Jenkins, A., Frenk, C. S., White, S. D. M., Colberg, J. M., 

Cole, S., Evrard, A. E., Couchman, H. M. P., & Yoshida, 

N. 2001, MNRAS, 321, 372 
Lacey, C, & Cole, S. 1993, MNRAS, 262, 627 
Lacey, C, & Cole, S. 1994, MNRAS, 271, 676 
Lahav, O., Lilje, P. B., Primack, J. R., & Rees, M. J. 1991, 

MNRAS, 251, 128 
Mo, H. J., & White, S. D. M. 1996, MNRAS, 282, 347 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 

490, 493 

Patiri, S. C, Betancort-Rijo, J., & Prada, F. 2006, MN- 
RAS, 368, 1132 

Peebles, P. J. E. 1980, Research supported by the National 
Science Foundation. Princeton, N.J., Princeton University 
Press. 

Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 

Reed, D. S., Bower, R., Frenk, C. S., Jenkins, A., & Theuns, 

T. 2007, MNRAS, 374, 2 
Sheth, R. K., & Lemson, G. 1999, MNRAS, 304, 767 
Sheth, R. K., & Tormen, G. 1999, MNRAS, 308, 119 
Sheth, R. K., Mo, H. J., & Tormen, G. 2001, MNRAS, 323, 

1 

Sheth, R. K., & Tormen, G. 2002, MNRAS, 329, 61 
Tormen, G. 1998, MNRAS, 297, 648 

Warren, M. S., Abazajian, K., Holz, D. E., & Teodoro, L. 

2006, ApJ, 646, 881 
Zentner, A. R. 2007, International Journal of Modern 

Physics D, 16, 763 



APPENDIX A: ANALYTIC EXPRESSIONS FOR 
THE MASS FRACTION OF DARK MATTER 
HALOES 

For some of the computations in this paper, we found that 
it is useful to have an analytical expression for the mass 
fraction F(m). We presen t here t hese expressions f or the 
case of IPress fc Schechteil ljl974h . ISheth fc Tormenl <(l999h 
and I Warren et alj (2006) mass functions. 

For the case of the PS mass function, we have the well 
known result of 

Fps(m) = erfcf 3, - ) (Al) 



\\/2a{m)) 

where erfc is the complementary error function, defined as 
erfc(x) = 2/0F/+ o °e- t2 dt. 

For the case of the ST mass function given in equation [§] 
we can integrate eq. [7] and obtain 



F ST (m) = A 



erfcf f 5c ) 



+ 



V^2p v 2 



~vA^f) (A2) 



where F(a, x) represents the incomplete Gamma function, 
defined as T(a, x) = f+°° e - t t a - 1 dt. 

Finally, for the case of the WA mass function given by equa- 
tion [10] we obtain 



-a/2 



(A3) 



where Ei(x) is the exponential integral function, defined as 



APPENDIX B: SOME USEFUL ANALYTIC 
FITS TO THE D FUNCTION 

We present here some useful fits to the D(q, Q, Q2) function, 
which was defined in eauation l24l All the numbers in this ap- 
pendix correspond to the transfer function and cosmological 
parameters presented at the end of Sec. 1 (see equation 1). 

In most of the cases considered in this paper, we are 
dealing with scales which are much smaller than the condi- 
tion (i.e. Q2 <S Q)- In that limit, the D(q, Q, Q2) function 



14 J. A. Rubino-Martin et al. 



is practically independent on Q2, and the following fitting 
function reproduces well the overall shape 



D(q,Q,Q 2 )=A(Q)e 



< ' 



(Bl) 



Note that in the limit q — > 0, we would expect this exponen- 
tial dependence to be exact. 

A numerical fit in the scale range 4 h _1 Mpc < 
Q < 13 h- 1 Mpc gives the following values: 



A(Q) = 1.264 + 0.167Q 8 - 0.0415Q|, 



(B2) 



lnB(Q) = -0.5632 + 0.1787(ln Q 8 ) - 0.0222(ln Q 8 ) 2 (B3) 

where Q 8 = Q/8 h" 1 Mpc. The typical error of this fit with 
respect to the exact computation is of the order of 5%. 

It is also interesting to provide a fit for the D function 
evaluated at the center of the condition, q = 0. In that case, 
we obtain 



D{q = 0, Q) = 1.262 + 0.1522Q 8 - 0.0392Q 8 



(B4) 



Finally, it is also useful to fit the function c(m) = 
D(Q, Q, Q), where m — m(Q). In this case, considering the 
mass range 10 8 h" 1 M < m < 5 x 10 15 h^ 1 M Q , we obtain 



i=4 

c(m) = S Ci (In m) 



(B5) 



with (c , ci, c 2 , c 3 , c 4 ) = (1.06004, -0.16521, 0.01024, -2.72x 
10 _4 ,2.96 x 10" 6 ). 



APPENDIX C: ANALYTIC FITS TO THE <5 NS 
FUNCTION 

In this section we provide analytic fitting formulae to the 
<5ns function, described in equation 1361 To a good approxi- 
mation, this function can be parameterised in terms of two 
functions, A(m, Q) and b(m, Q), which are defined in equa- 
tion [38] 

All the following fits have been obtained with the ST 
UMF and using 8 e g = 1.25. We have considered values of Q 
within the range between 5 and 13 h _1 Mpc. The mass in- 
terval for all fits is taken to be 10 9 h _1 M Q < m < m(Q)/3Q. 
However, we find that the functional form proposed in equa- 
tion [38] only provides a reasonable fit (i.e. with errors of few 
percent) to the data for masses above m/m* > 3 x 10 -3 
(where m* is defined as <r(m*) = 1). For mass values outside 
that range, it is necessary to do the numerical integration 
in Eq. [36] Finally, all these fits have been obtained for val- 
ues of the linear density within the range —4.5 < 5\ < — 1. 
These are the typical values we are interested in. We note 
that beyond this range, and in particular, for values of 5\ 
close to zero, these fits may give inaccurate results. 

The proposed fitting formula for A(m) is simply a 
quadratic fit in him, i.e. 



A{m,Q) =^a 1 (Q)(lnm') 1 

where m' is the dimensionless mass given by 
, m 



3.51 x 10 11 h" 1 M Q 
and the coefficients are linear functions of Q, 



(CI) 



(C2) 



a,i(Q) 



if Q 

Hsh-i 



Mpc 



(C3) 



The results we obtain are (ao,Oo) = (1.577,-0.298), 
(a1,al) = (-0.0557,-0.0447), and {a 2 ,a\) 
(-0.00565,-0.0018). 

The best fit for b(m) is obtained with the formula pro- 
posed in lPatiri et~ai1 l|2006l ). 



b(m,Q)^b 1 (Q) + b 2 (Q) (m') bM) (C4) 
but the coefficients are now linear functions of Q, 



(C5) 



Using these expressions to fit the numerical results, we find 
(61,61) = (-0.0025,0.00146), (62,63) = (0.121,-0.0156), 
and (63,63) = (0.335,0.019). 



APPENDIX D: ANALYTICAL FIT TO THE 
CONDITIONAL MASS FUNCTION FOR ANY 
COSMOLOGY 

Following the approach presented in this paper, we give here 
a highly accurate analytical fit to the CMF, as a function 
of the redshift and the cosmology. The (accumulated) con- 
ditional mass function N c (m), as a function of the linear 
density contrast 5i, within a spherical conditioning region 
with Lagragian radius Q is given, for any value of the red- 
shift z, and the cosmology (cr 8 and fi m ), by 



N cE (m, z\Q, 8 l ;a s ,F)= Ae~ bd i (1 + S)N(m, z), 



(Dl) 



where N(m, z) is the (unconditional) accumulated mass 
function, S is the actual fractional mass density fluctuation, 
and the coeffitients A and b are functions of m, Q, z, u 8 
and r given by 



A(m, Q; z, as, T) 



b(m,Q; z,as,T) 



asb(z) 
0.9 



o s b(z) 
0.9 



0.71 + 0. 08m' 



A(m,Q) 



6(m, Q) 



(D2) 



(D3) 



and where b(z) is the growth factor of linear density per- 
turbations normalised to one at 2 = 0, and A(m, Q) and 
b(m,Q) are given by equations ICll and lC4l respectively, but 
using this new mass definition 



(D4) 



3.51 x 10 11 h- 1 M, 




It is found that the V dependence of both A and b 
is negligible, but the dependence on <r 8 (and therefore on 
redshift z) is rather strong, particularly in the case of b. 
These two scalings presented in equations ID2I and ID3I have 
been obtained as a fit to the numerical results for values of 
the linear density in the range —4.5 < <5i(0.9/<7g) < —1, as in 
the previous appendix. We have checked that changing the 
range of integration may produce changes on the exponent of 
these scalings up to the order of 10 per cent. The exponent 
in both cases has been fitted in the mass range m > 3 x 
10- 4 m*. 



Finally, it must be noted that the range of values of as 
and F used to obtain this fit is roughly a factor two around 
the reference model (i.e. as = 0.9 and F = 0.21). Thus,it can 
not be safely extrapolated beyond z w 2. Furthermore, for 
la rge redshift values (z > 5), the considera tions contained 
in iBetancort-Riio fc Montero-Dortal |2006b) must be taken 
into account. 
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